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We apply Euclidean time methods to phenomenological Dyson-Schwinger 
models of hadrons. By performing a Fourier transform of the momentum 
space correlation function to Euclidean time and by taking the large Euclidean 
time limit, we project onto the lightest on-mass-shell hadron for given quan- 
tum numbers. The procedure, which actually resembles lattice gauge theory 
methods, allows the extraction of moments of structure functions, moments 
of light-cone wave functions and form factors without ad hoc extrapolations 
to the on-mass-shell points. We demonstrate the practicality of the procedure 
with the example of the pion form factor. 

I. INTRODUCTION 

The planned experimental program at the Thomas Jefferson National Accelerator Facil- 
ity (TJNAF) will subject the electromagnetic structure of hadrons and nuclei to detailed 
scrutiny. The calculation of hadron form factors is therefore of fundamental importance to 
the subsequent interpretation of the obtained experimental results. However the nonpertur- 
bative description of hadron structure in the Minkowski metric, which has only recently been 
pursued in detail , is extremely difficult due to the direct confrontation with singularities 
and the indefinite norm. Alternatively, Euclidean space is characterized by a positive definite 
norm, i.e., p 2 = p\ + p\ + p\ + pi > 0, and it has long been known that the Euclidean for- 
mulation is therefore advantageous in the description of nonperturbative processes through 



the Dyson-Schwinger(DSE) and Bethe-Salpeter(BSE) equations. 

These components are frequently assembled into diagrams such as that in Fig.l for 
the calculation of form factors (0,0. The problem with the Euclidean formulation of such 
processes is that for physical particles the external momenta must satisfy the mass-shell 
condition P 2 = —M\ and K 2 = — Mf , thus forcing the return to Minkowski space. One is 
then faced with the problem of complex momenta flowing through the loop in the diagram 
of Fig.l. This in turn requires solving the DSEs in the complex plane. Although some 
progress has been made in this regard @, it is an extremely difficult problem and most 
calculations until now have employed entire-function fits on the real axis. Although some 
success has been obtained with this approach ||, there are uncertainties associated with the 
extrapolation of these functions into the complex plane. 




FIG. 1. The electromagnetic vertex for a composite pion at tree level. The quark Green's 
functions, quark-photon vertex, and the pion Bethe-Salpeter amplitudes, are dressed in a consistent, 
gauge-invar iance manner. 

Here we offer a fresh approach which avoids these uncertainties by obtaining the mass- 
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shell conditions implicitly through an application of the Cauchy integral formula]] This 
approach does not circumvent the need to make assumptions about the analytic structure 
of the components in the diagram of Fig.l, but rather relies on the selection of relevant 
singularities in much the same way as is accomplished by Euclidean-space lattice calculations. 
We consider the pion form factor as a prototype for the purposes of illustration. 

The exact calculation of the pion form factor via lattice techniques has been studied 
previously and proceeds as follows: one first considers the Euclidean three point cor- 
relation function 

C,(t x ,t y ) = J d 3 xd 3 y(0\T [jl(y)J,(0)Ux)} |0)e" ii? 'V^, (1) 

where J n (x) = u(x)i^d(x) is the interpolating field for the pion and J M (0) = |-u(0)7 At n(0) — 
|J(0)7 M (i(0) is the vector current. t x and t y are the time components of x and y respectively. 
Upon inserting a complete set of states between each pair of operators in Eq.(|l|) and taking 
the limit t x — > — oo and t y — > oo, so that only the lightest state contributes significantly, one 
finds 

, lim C^ty) = - 6 " (014(0)1^^)^,^1^(0)1^^)^,^1^(0)10), (2) 



where \n,P) represents an on-mass-shell (i.e. Ep = \J M% + P 2 ) pion state with three- 
momentum P. Note that Eq.@ is exact. Thus despite the Euclidean formulation (as long 
as one does not employ any further approximations such as quenching, finite quark masses, 
finite lattice spacing, finite volume) it allows an exact computation of the on-shell pion 
electromagnetic form factor 



(P„ + Kfj)F^ UP - Ky) = (tt, K\J,(0)\n, P), (3) 



1 A similar method has been applied in Ref. || to 2-point functions, where it has been used as 
a means to obtain the physical mass. However, in this paper, we would like to demonstrate that 
Euclidean time projection methods are much more powerful and can be successfully applied to 
higher point functions as well. 
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where P = (Ep,P) and K = (Eg,K). In practice 



Eq. QT|) is usually evaluated in the 



quenched approximation, yielding 




dUe~ s{u) J d 3 xd 3 ye 




tr 75 M- 1 ([/)(x,0)7 M M- 1 (f/)(0, ? /)75M- 1 ([/)( 2 /,a;) 



(4) 



where M~ 1 (U)(y, x) is the quark propagator for a given configuration of link fields U. Note 
that the group integration / dU implies that the quark lines and the vertices are dressed 
with arbitrary gluon lines. 

The above example (plus many other examples in Euclidean lattice gauge theory) shows 
that, despite the fact that the calculations are performed in a Euclidean metric, it is still 
possible to extract on-mass-shell matrix elements, without having to make ad hoc assump- 
tions about the analytic continuation back to a Minkowskian metric. The basic idea is that 
large Euclidean time merely acts as a filter to project out the ground state of the Hamilto- 
nian. In this paper, we present a procedure that uses similar techniques in the context of 
Euclidean Dyson-Schwinger calculations, which also there allows an unambiguous extraction 
of the on-mass-shell matrix elements. 



An alternative to performing the exact (or quenched) calculation is to employ two ap- 
proximations restricting the type of allowed gluon dressing. The first is to consider only 
dressing by the gluon two-point function, and the second is to consider only ladder /rainbow 
dressing .0 These approximations are in fact implied by the diagram of Fig.l, and allow 
the separate calculation of the components, i.e., quark propagators and vertex functions. 

2 Note that ladder /rainbow dressing implicitly leads to a restriction on allowed quark lines that 
resembles the quenched approximation in lattice QCD. 

4 



II. THE DYSON-SCHWINGER APPROACH TO FORM FACTOR 



CALCULATIONS 



These approximations further comprise an electromagnetic gauge-invariant description 
Within this approximation the diagram of Fig.l is given by 



Al(K, P) 



d A k 



tr 



A 5 (P,fc)G(fc+^)A A , (K-P,k+j 



(5) 



where A 5 is the solution of the inhomogeneous pseudoscalar BSE (quark-pseudoscalar ver- 
tex), A M is the solution of the inhomogeneous vector BSE (quark-photon vertex), and G is 
the quark propagator. 

We will use the notation A to represent the solutions to the inhomogeneous BSE, and T 
to represent the solutions to the homogeneous BSE. In the vicinity of a pole, they are related 

by i 



A(P,k) 



M 2 JZ(-M 2 ) 



r(p,k), 



(6) 



P 2 + M 2 

where M is the pole mass, Z(-M 2 ) is the wave function renormalization, and T(P,k) 
therefore obeys the standard Bethe-Salpeter bound state normalization ||. 

The vertex function for a photon coupling to an on-shell pion is then given by 



r dk 
K(K,P) = / t — -r tr 



(27 



where in the vicinity of the pion pole 



T n (P, k)G(k + |) A„ (K-P,k+^- 
xG (k - - + k)t v (K,k - ^ + — ) G (k - - 



(7) 



MtZ(-Ml) 



(P 2 + M 2 )(K 2 + M 2 ) 



K(K,P). 



(8) 



The quantity T*(K,P) has the decomposition 



r;(K, P) = {P, + K,) F n (P 2 , K 2 ,(K — Pf 

+ (K^ - P,){P 2 - K 2 )H n (P 2 , K 2 , (K - P) 2 ) (9) 

and is commonly used to extract the pion form factor, F n , at the on-shell point P 2 = K 2 = 
-Ml. 



Rather than explicitly enforcing the mass-shell condition as has been explored previously 
||, here we instead take an approach that is similar to that taken in the Euclidean-space 
lattice calculation. The quantity in the DSE approach that is analogous to the three-point 
correlation function, C^{t x ,t y ) in (|J), is 

Cl(t x ,t y ) = -j~J^± e-^e^Al(K.P). (10) 

The form factor is then obtained implicitly from the expression 



lim C'{t x ,t y 



M*JZ{-Ml) 



The correctness of Eq.([0]) is readily verified by substituting Eqs.(|]) and (|) into QlCf), along 
with the implicit assumption that there are no poles in the immediate vicinity of the ground- 
state pion pole of interest. The validity of this assumption is of course model dependent. 
Eq.flnp can be directly compared with the exact result in Eq.(|2|). 

III. TOY-MODEL CALCULATION 



We now proceed to evaluate Eq. (|Tl|) for a simple model which allows direct comparison 
with the exact result. That is, we make the following choices 

A^P,k) = -i llx (12) 

G(k) = -. — r 

27 • k + m q 

where m q is taken to be a constituent quark mass. Eq.(|12"D implies in particular that = 
i / y 5 . The exact result can therefore be obtained by directly fixing the mass-shell condition in 
the evaluation of Eq.(|7]) and Eq.(|9|), and is shown by the solid line in Fig. 2. The constituent 
quark mass, m q , is taken to be 300 MeV, and a cutoff of 1 GeV is applied on the relative 
momentum at the quark-photon vertex to regulate the integration. The external momenta 
are constrained such that P 2 = K 2 = 0. 



For comparison, we evaluate Eq.(|TT|) by substituting the toy model ( |12"D into Eq.([5D 
for the vertex A^. The results are shown by the series of dashed curves in Fig. 2 for time 
separations At = t y — t x ~ 2fm. The Fourier transform is performed with 4096 equally- 
spaced points with a grid spacing of 0.02 GeV. Note that the agreement with the exact 
result diminishes at larger momentum transfers Q 2 , as is expected. The agreement at large 
momentum is improved by increasing the number of Fourier transform points to include 
larger momentum components in the integration. 

There is an additional effect due to the presence of poles in the quark propagators. 
This is revealed in Fig. 3 where for example a constituent quark mass of 1 GeV is used 
with 1024 Fourier transform points at a grid spacing of 0.03 GeV. With this large quark 
mass the associated singularities are strongly damped and therefore do not contaminate the 
extraction of the pion pole. The comparison with the exact result is now quite good out 
to reasonably large momenta. The increasing discrepancy with increasing time separation 
at large momenta is due to the highly oscillatory nature of the Fourier transform, and is 
improved with a finer grid spacing. It is expected that the use of quark propagators from 
DSE studies || that reflect the confining nature of QCD will also produce accurate results 
at higher momenta due to the absence of such poles. 
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FIG. 2. The exact and extracted "pion" form factors are compared for a constituent quark mass 
m q = 300 MeV. 
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FIG. 3. The exact and extracted "pion" form factors are compared for a constituent quark mass 
m q = 1.0 GeV. 
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IV. SUMMARY AND OUTLOOK 



We have introduced Euclidean time projection methods as a means to project Euclidean 
metric 3-point functions onto the physical mass shell. The method, which is based on 
Cauchy's theorem and the existence of a Lehmann representation, avoids having to make 
ad hoc assumptions about the analytic behavior of propagators and vertex functions, and 
thus allows the direct use of Euclidean DSEs for the calculation of form factors. We have 
demonstrated the feasibility of the method by studying the pion form factor for a toy model 
where the exact form factor is known. The approach resembles methods that have been used 
in the context of lattice gauge theory The idea is to first Fourier transform n-point 

functions from (Euclidean) energy to (Euclidean) time. The Euclidean time difference be- 
tween hadron interpolating fields and currents that are used to "probe" the hadron structure 
is then taken to be large. This ensures that the hadron which is being probed is actually 
the ground state with the quantum numbers of the interpolating field. 

Even though we have studied the pion form factor as a showcase example, applications 
of Euclidean time projection methods should be applicable to a much wider class of physical 
observables. Starting with 2-point functions, one can obtain masses || as well as light-cone 
moments of hadron wave functions J7| on the mass shell. Examples for 3-point functions that 
we plan to calculate include the elastic form factor, moments of parton distribution functions 
(e.g. momentum or spin fraction carried by the quarks in a given hadron), and the 7*7r — > 7 
transition form factor. The most simple observables related to 4-point functions that one 
might consider calculating using Euclidean time projection are hadron polarizabilities and 
(virtual) Compton scattering processes. 

The main limitations of the method are first of all the restriction to the lightest hadrons 
for given quantum numbers: though theoretically possible, it is in practice rather difficult 
to use the method to obtain observables that involve excited states, since the whole trick 
is to use large Euclidean times to suppress excited states (a limitation that is very familiar 
to the lattice gauge community). Another limitation concerns very large momentum trans- 
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fers, where rapid oscillations of the integrand require a large number of integrations points. 
However, even with rather moderate computing power, this still allows access to most of the 
kinematic range that is available at TJNAF. 
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